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Abstract: We use large-scale numerical simulations to study the formation and evolu- 
tion of non-topological defects in a generalized electroweak phase transition described by 
the Glashow-Salam- Weinberg model without fermions. Such defects include dumbbells, 
comprising a pair of monopoles joined by a segment of electroweak string. These exhibit 
complex dynamics, with some shrinking under the string tension and others growing due to 
the monopole-antimonopole attractions between near neighbours. We estimate the range 
of parameters where the network of dumbbells persists, and show that this region is nar- 
rower than the region within which infinite straight electroweak strings are perturbatively 
stable. 
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1. Introduction 

The formation of topological defects in any phase transition depends on the existence of 
a non-trivial, low-order homotopy group of the broken-symmetry vacuum manifold |, §. 
Even in the absence of this, dynamically stable non-topological defects are sometimes still 
possible, but are usually assumed to be too weakly stable to lead to a lasting network. A 
particularly interesting case is non-topological strings where at least one example is known 
(semilocal strings [Q, Q ) in which a space-spanning network of strings can develop from the 
growth and joining of short string segments formed at the phase transition ||. Semilocal 
strings are a special case of electroweak strings |], 0, ||]. 

The stability of electroweak strings has been analyzed in some detail in the Glashow- 
Salam- Weinberg (GSW) model. In the absence of fermions the theory is determined, up to 
scalings, by only two parameters: the weak mixing angle an d which is the square of 
the ratio of the Higgs and Z-boson masses (the observed value for the actual electroweak 
model is sin 2 6>w — 0.23, and the precise value of (3 has not yet been determined but is 
almost certainly greater than one). 

In the case of interest here, the topology of the vacuum manifold is the three-sphere 
S 3 , which does not support persistent topological defects in 3+1 dimensional spacetime, 
but there are grounds for expecting non-topological defects to form. In the limit #w = 
7r/2 we recover the so-called semilocal model, whose S 3 vacuum manifold is known to 
support stable non-topological strings from both analytical [||, Q and numerical Q work. 
In particular, a straight infinite semilocal string can be seen as a Nielsen-Olesen (NO) U(l) 
vortex @, embedded in a higher group SU(2) g i b a i x U(l)i oca i. Note, however, that in the 
semilocal string case, we can define a quantity (the winding number) which, although not 
a topological invariant in the usual sense, is topologically conserved. 



The embedding of NO vortices into the full electroweak symmetry SU(2) l x U(1) y leads 
to Z-strings and VT-strings (see ref. ||). Although VK-strings are expected to be unstable, 
analysis of infinite axially-symmetric Z-strings has shown the existence of a parameter 
regime where they are perturbatively stable |7|, [T0| . But these Z-strings are genuinely 
non-top ological, there is no quantity which is topologically conserved. 

In this paper we will consider the evolution and persistence of these genuine non- 
topological defects in a generalized GSW model, spanning all values of /3 and 6*w, using 
numerical methods. Isolated, infinite, axially-symmetric Z-string configurations will not 
be formed in a realistic system, but configurations consisting of monopole-antimonopole 
pairs joined by Z-string segments, named dumbbells by Nambu ||, are certainly possible. 
Isolated dumbbells are expected to collapse under the string tension, at least in the absence 
of rotation, jittering or magnetic fields [||, llj 12]. Of particular interest is the question of 



whether there is a region of parameter space in which the density of dumbbells is sufficiently 
high so that they are able to generate a persistent network of strings by building up longer 
segments from the merging of shorter ones due to interaction of neighbouring monopoles, 
as has been found in the semilocal case ||. This region of parameter space will lie far 
from the measured values for #\y and 0, but the mere existence of this region is striking. It 
shows that in models close to real physical ones, completely non-top ological defect networks 
can form and persist, so in extensions of the GSW model (to higher symmetry groups or 



extra fields), or models with a background magnetic field [12], or even models in which 
topological (semilocal) defects exist in some limit, non-topological defects cannot be ruled 
out immediately. 

Moreover, the study of the dynamics of dumbbells can be of interest in computing 
other early universe features, such as the primordial magnetic field helicity |l5|, O ], 

2. The model 

The bosonic sector of the GSW electroweak model describes an SU(2)^ x U(1) Y invariant 
theory with a scalar field <& in the fundamental representation of SU(2)i, with lagrangian 



C = \D^\ 2 - ^W^W a ^ - - A (W - ^] . (2.1) 



The covariant derivative is given by 

D^d.-^r^-^Y,, a = 1,2, 3, (2.2) 

where $ is a complex doublet, r a are the Pauli matrices, W£ is a SU(2) gauge field and 
Y^ is a U(l) gauge field. The field strengths associated with these gauge fields are 

W%, = W - d v W« + g w e abc W b ^ ; 

Y^ = d^Y u - d u Y^ , (2.3) 



respectively, and there is no distinction between upper and lower group indices (e 123 = 1). 



When the scalar field acquires a non-zero vacuum expectation value the symmetry 
breaks from SU(2)i x U(1) Y to U(l) e . m . , leaving a massive scalar field (m# = y/2Xrj), a 



massless neutral photon (Au), a massive neutral Z-boson (Z ( 



A" 



g z r)/2 = l v 1 , where 



= V 9y + ffw); an d two massive charged VF-bosons (Wjf , mvK = 9w r ?/2)- 
We make the following rescaling 



v 7 ^ 



V2 
x 

9zV 



9w w* - 9 -^w; 



(2.4) 



to choose l v as the unit of length, r] as the unit of energy, and the Z-charge of the scalar 
field (g z ) as the unit of charge (up to numerical factors). This brings the classical field 
equations to the form 



2A / \ 
D^D^ + -j ($$ f - lj $ = : 



d v W' J *' a + e abc W^W^ vc = ^ 



cos #w 



= ^sin 2 # w $ t r> M $- (-D M $)* $ 
where the weak mixing angle is given by tan#vv = 9y/9w, an d now 



(2.5) 



W a 
D„ 



2 M 



-Y„. 



(2.6) 



Typically the Z- and A-fields are expressed in the unitary gauge <& T = (0,1), but this 
choice is not well suited to work with defects. Instead, when there are points in space-time 
with |<3?| ^ 1, it is customary to use a more general definition of these fields which depends 
on the Higgs field configuration at each point Q, namely, 



= cose w n a (x)W l 



A, 



p sin W ; 
sin 6>w n a (x) + cos 6\v^ , 



where 



n a (x) 



$t( x ) T a$( x ) 
&(x)<S>(x) 



(2.7) 



(2.8) 



is a unit vector by virtue of the Fierz identity ^ a (<J>V a< I>) = ($^$) . There are also 
several possible definitions for the field strengths (see, e.g. ref. |l3| for a discussion of this 
point); in our simulations we will use 



Z^ = cos% n a (x) W^ v - sin% Y^ v ; 
A^ u = sin% n a (x) W" + cos% Y, 



■ [IV 



(2.9) 



Note that eqs. ( |2.7[) to (|2.9| ) reduce to the usual definitions away from the defect cores. We 
work in flat space and in the temporal gauge (Wq = Yq = for a = 1, 2, 3) so Dq<& = do<&. 
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With this gauge choice, Gauss's Law becomes 

i 
2 



djidoYj) = ^sin 2 w [* f Sb* - (do$) f 3> 



-dj{d Wf) - e abc W}d W? = ^cos 2 9 w |$V a 0b$ - (<9 $) T r a $j , (2.10) 

(with j = 1,2,3) which is then used to test the stability of the code. 

As noted above, the sin 2 #w = 1 case reduces to the semilocal case, where absolutely 
stable defects are known to exist. Setting one of the Higgs fields to zero further reduces 
this to the abelian Higgs model. We are able to use our knowledge of these systems, as well 
as the behaviour of artificially constructed infinite axially-symmetric Z-strings at various 
points in parameter space, to check the validity of our simulation code. 

3. Numerical simulations 



The equations of motion eq. (2.5) are discretized using a naive staggered leapfrog method, 
where both the scalar and gauge fields are associated with lattice points. This procedure 
is the easiest to relate to previous work on semilocal defects ||, and we expect to have 
a fairly good picture at the scales we are interested in. A comparison to a link-variable 
discretization of eq. ( f2~5| ) in the semilocal case (#w = 7r /2) shows that the changes 
from using a lattice implementation instead of this naive one are well within the other 
uncertainties; see the appendix for a full analysis. The simulations are performed on a 
periodic cubic lattice whose time step is 0.2 times the spatial step At = 0.2Ax (c = 1), with 
ad hoc numerical viscosity terms added to each equation (7$, jY and "fW a respectively) 
to reduce the system's relaxation time. The expansion rate in an expanding universe 
would play this role, albeit as a time-dependent factor ^(t) typically scaling as 1/t. Several 
different values of 7 were tested, and generated very similar behaviour; throughout this 
work we use 7 = 0.5. 

Two different strategies for setting the initial configurations were considered: 

(a) Set all initial field velocities to zero and throw down random scalar field phases at 
every lattice point; average the field at each point with its 6 nearest neighbours 
and normalize the fields to the vacuum value (|$| 2 = 1 as r\ has been rescaled out) 
iteratively (50 times), to get a smoother configuration; using these scalar field values, 
choose the initial values of the gauge fields to be 

Y„ = 

Wl = 2 (ViV^4 - faVjih + tfeVj^j - VaVjVa) 

Wl = 2 (VsV^i - V'iVjVs + faVjih - faVjfa) 

W* = 2 (^iV^2 - VaVjVi + V^V^s - V'sV^) (3.1) 

where <1> T = (ipi + i ip2 > ^3 +*V'4)j pseudo-minimizing the energy as in ref. |L4|]. Note, 
however, that in the present case, the system does have enough gauge fields to cancel 
completelly the gradient energy. In fact, the configuration ( |3.1|) sets to zero the gra- 
dient energy, the potential energy and the field-strength energy corresponding to the 
U(l) gauge field, leaving the SU(2) field-strength energy as the only non-zero term. 



(b) Set all fields to zero, and also gauge field velocities to zero ($ = W£ = = = 
= 0); provide some initial (smoothed) random velocities to the scalar field (<I> 7^ 0) 
following the general procedure as above. This initial configuration would be closer 
to that appropriate to a bubble nucleation scheme. 

The overall results in our simulations did not show qualitatively different behaviour 
between the two cases. Indeed, in the first case the gauge field energy pseudo-minimization 
was ineffective enough that the scalar field typically began by climbing up the potential, 
restoring the symmetric phase, and rolled down again later. Previous experience with 
semilocal strings also shows that the results are fairly insensitive to the specific way initial 
conditions are implemented. For all these reasons, we chose the second initial configuration 
(6) for our investigations. 

Interpreting simulations of electroweak string networks is more complicated than in the 
U(l) cosmic string case. As in the semilocal case, electroweak strings are non-topological, 
and in this case there is no well-defined winding number, which makes identification of 



strings a more difficult task. We follow the strategy proposed in ref. [14] in order to 
study string formation: the system is evolved forward in time and we compute a set of 
gauge-invariant quantities at each time step, namely the Z- and the ^4-field strengths 



given by eq. (|2.9|) and the modulus of the scalar field (|3>|). String formation is then 



observed by visualizing isosurfaces in the Z-field strength (W ^ZijZ l i) and in the scalar 
field modulus. 

After the rescaling eq. (|2.4| ), it becomes clear that the only free parameters in our 
model are f3 (/3 = m 2 H /m1 = 8\/g%) and #w It is well known that with Qyj = it/2 
and (3 < 1 we have semilocal strings that are stable 0, and numerical simulations show 
segment formation and linkage [||. Moreover, for sin 2 ^ ^ 1, /? < 1, there is a regime, 
albeit a rather narrow one, where infinitely long strings are perturbatively stable |7|, 10]. 



Bearing these results in mind, the parameter space investigated in our simulations was 
0.9 < sin 2 w < 1 and 0.05 < (3 < 1.5. 

Simulation and visualization were performed on 64 3 , 128 3 and 256 3 lattices using the 
Cray T3E at NERSC, and high-performance computing facilities at University of Sussex 
and University of Groningen. Our quoted results all come from 256 3 simulations. 



4. Results 



To test our code, we began by reproducing known results for cosmic and semilocal string 
networks, and for infinitely long, axially-symmetric Z-strings (made possible by our peri- 
odic boundary conditions as long as strings are simulated in pairs to keep the net flux equal 
to zero) in both the stable and unstable regimes, in particular verifying that the Z-strings 
disappeared in the unstable regime. 

Having checked the code, we then ran it for values of the two free parameters (9w,P) 
throughout the regime of interest, using the same initial conditions for each parameter pair. 
We chose to carry out one large simulation at each point in parameter space rather than 
many smaller ones, which gives improved dynamical range. As expected, after an initial 



Figure 1: Isosurfaces of the Z- and ^4-magnetic field strengths for two different simulations, shown 
as the light and dark colour respectively. The left panels show (3 = 0.1, sin 2 6>w = 0.994 (persistent 
regime) and the right ones /3 — 0.5, sin 2 #\y = 0.995 (non-persistent regime). The top row is at an 
early stage of the simulation t = 50, while the lower is at the end t — 200. Note that in the first 
case there remain some long strings at the end of the simulation and connection can still occur. In 
the second all the defects are about to disappear. 

transient, the typical configuration observed is a dumbbell — a segment of electroweak 
string joining a monopole/antimonopole pair. By continuity, for parameter values suffi- 
ciently close to the stable semilocal case we expect the monopole interactions to lead to 
the joining of strings to form longer segments. 

The simulations show that some short segments of Z-string do join as expected. The 
joining rate is, however, lower than in the semilocal case, and decreases both as sin 2 #\y is 
decreased and as (3 is increased. Thus, some segments which would eventually join in the 



-6- 



semilocal case are seen instead to collapse in the electroweak case due to string tension. 
This is not surprising: in the semilocal case, we have global monopoles at the string ends, 
which have divergent scalar gradient energy and are more efficient at finding neighbouring 
monopoles. In the electroweak case, the monopoles at the string ends are proper magnetic 
monopoles, and the scalar gradients are cancelled much more efficiently by gauge fields. As 
sin 2 #w — ► 1 the cores become larger, and eventually overlap, making segments join. But as 
we move away from the semilocal case, the cores become smaller and the joining becomes 
less important. 

Figure |l| shows timeslices of two typical simulations, with the different colours corre- 
sponding to the A- and Z-field strengths. 1 

The Z-field has a string-like form, whereas the A-field at the string ends is a spherical 
shell, corresponding to spherical magnetic monopoles. The ^4-field morphology can also be 
tube-like, denoting interaction between monopoles, illustrating the complexity of the overall 
dynamics. The interplay between string tension and monopole-antimonopole attraction 
causes some strings to shrink until they disappear, and others join to form longer strings. 

To compare the number of defects in the semilocal and electroweak cases, we plot the 
number of lattice sites with Z-magnetic field strength (^J ^ZijZw) greater than 25% of the 
maximum field strength found in the core of a NO string with the same value of (5 (see 
figure ||). The number of lattice sites, and hence length of string, decreases more rapidly 
in the electroweak case as either sin 2 #-w decreases or (3 grows. 

At early times the configurations in different simulations are very similar. During the 
symmetry-breaking transient (the very first time steps) no defects exist. As the scalar field 
takes on a non-zero value a large number of very small string segments emerge, then the 
magnetic monopoles become visible at the ends of the segments due to the gathering of 
A-magnetic flux there. The upper panels of figure |l] show the Z- and ^4-magnetic field 
strength of two different simulations at time t = 50. Similarities in these first time steps, 
together with similarities using different initial conditions (see above) , show that the initial 
configurations (i.e. the way the symmetry breaking is implemented) are not as important 
as the subsequent interaction between the scalar and gauge fields. 

As the system evolves, we see that in only one of the two simulations do the small 
segments grow and connect to their neighbours, allowing the defect network to persist. In 
the lower panels the configuration at time t = 200 can be seen. In one case there are long 
strings and connections are still happening, whereas in the other almost all of the dumbbells 
have annihilated. In particular, in the first case the final configuration contains strings 
which are much longer at the end of the simulation than any present in the early stages. 

We wish to determine the parameter space for which the string network persists. The 
stability of the corresponding infinite Z-strings is a necessary condition for persistence, but 
not sufficient since the physically realizable collection of finite initial string segments may 
not link up. We therefore expect the persistence region of parameter space to lie entirely 
within the stability region. 

1 Further colour images and movies can be found at 
http : //www .nersc . gov/ ~borr ill/defects/ electroweak.html 
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Figure 2: The number of lattice points with Z-field strength (y ^Z^Z^) bigger than 25% of the 
maximum (core) value in the NO case, calculated in a series of 256 3 simulations, for j3 = 0.3. The 
different lines correspond to different choices of the weak mixing angle #\y (values shown are for 
sm 2 6\y)- Note that in our criterion for persistence (see text), for the chosen f3, the defects are 
considered to live long enough for sin 2 #w ^ 0.995. Different initial conditions for the simulations 
give similar behaviour. 



To estimate the persistence of Z-strings we need to define a criterion to decide when 
defects are lasting long enough. There is clearly some degree of arbitrariness in how this 
criterion is set, and we have been guided by inspecting visually the evolution in different 
parameter regimes. In order to construct the persistence region in parameter space we 
simulated the system in 256 3 boxes with periodic boundary conditions, and considered 
persistence to have occurred if at time t = 200 there are more than 1000 lattice sites with 
a Z-magnetic flux greater than 25% of the maximum of the NO simulation for the same 
(3. For example, for f3 = 0.3, we can see in figure || that persistence defined this way is 
exhibited only for sin 2 #w > 0.995. From a suite of simulations using this criterion we obtain 
the 'persistence limit' shown in figure [3|. As anticipated, the persistence region covers only 
a subset of the stability region, with persistence only for values of sin 2 #w extremely close 
to one. Another possible criterion for persistence which can be easily automated is the 
presence at late times of strings that are reasonably long compared to their width. 2 At 
time t = 200, the condition that there be strings at least five times their width leads to 
the second persistence line drawn in figure ||. 

2 To automate the calculation of the length of the individual strings, we compute the volume of each 
segment of string by counting the number of connected points whose Z-magnetic field is at least 25% of 
the maximum in a NO string at the same value of (3. Then, the volume is divided by the cross-sectional 
area of the NO string. Finally, to obtain the length-to-width ratio, we divide it by the diameter of the 
corresponding NO string. The radius of the NO string is always taken to be the distance between the center 
and the point where the field strength is 25% of the maximum. 
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Figure 3: The solid line is the semi-analytical curve marking the stability transition for infinite, 



axially-symmetric Z-strings 10 , and the points mark the edge of the persistence region obtained 
in our simulations, using both criteria described in the text. As expected, the persistence limit 
lies entirely within the stability zone. The squares are obtained considering that the defects are 
persistent if at t — 200 there are more that 1000 points with a magnetic field strength higher than 
25% of the maximum (see text); the triangles are obtained by considering persistent those defects 
which at t = 200 are at least 5 times longer than their width. 



5. Conclusions 



The numerical simulations described here show that a significant non-topological defect 
network can form in a generalized GSW model for an electroweak phase transition. The 
dynamics of such a network are extremely complicated, driven by string segment Unkings 
and by isolated strings shrinking, and the details are highly sensitive to the two model 
parameters #\y and (3. Our principal result is that in some regions of parameter space a 
persistent network of genuine non-topological defects can form. Though the actual version 
of electroweak theory in our Universe lies outside this parameter regime, these results show 
that, in models where topological (or semilocal) defects are possible in some limit, it is 
possible to get a network of non-topological defects close to that limit. Previous works in 
the literature [jj], 1C , ^] show that for parameters close to those permitting topological (or 
semilocal) defects, there is a regime where non-topological defects are stable. Our work 
confirms this result, and shows that, although in a narrower region, a sufficiently persistent 
network of defects can form in such a phase transition. 

Note that the generation of the string segments is intrinsically dynamical and cannot 
be studied using initial condition arguments; in particular our results are compatible with 



those of [18]. Here we are considering the time evolution of a network and not just the 
initial configuration. In fact, the first timesteps in such simulations correspond more to a 
numerical transient, in which physically reasonable initial conditions are established, than 
an actual phase transition. Only after this initial transient can the evolution of the network 



Figure 4: Plot of the magnetic field strength for the semilocal case at t = 40 (top figures) and 
at t — 56 (bottom figures), using the same initial conditions, but evolving the system using a naive 
discretization (left) and a lattice link variable discretization (right), as explained in the appendix. 
Both simulation schemes are in very good agreement, with slightly longer strings in the lattice 
variable discretization scheme. 



be trusted, and this is what determines whether the defects persist or decay. Given the 
primary role played by the gauge fields, it would be interesting to know if our conclusions 
generalize to non-top ological string defects in other models such as the two-Higgs standard 
model MM. 
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A. Discretization methods 



In this work, we discretized the equations of motion eq. ( |2.5| ) replacing all the scalar and 
gauge fields by their values at the lattice points. For instance, a U(l) covariant derivative 
was substituted by 

n il \ 4>(x + Xj) - (j)(x - Xj) igy^ rr \ i / \ ,k n 

D j( j){x) = ^ ^ - ^-Yi(x)4>(x) (A.l) 

where L is the lattice spacing. 

It is well known that this discretization scheme does not respect gauge invariance. 
Nevertheless, we are dealing with well-defined classical equations of motion in a particular 
gauge, and the breaking of gauge invariance should be irrelevant in this context. Indeed, 
we make use of the fact that Gauss's Law is not automatically satisfied, and by monitoring 
it during the evolution we can check that our discretization is reasonably accurate. This 
discretization was chosen to facilitate comparison with earlier work [jlj| ]. 

There is an alternative discretization method widely used in the literature which pro- 
tects gauge invariance and recovers the original equations of motion in the limit where the 



lattice spacing goes to zero [17]. This method uses lattice link variables, i.e., the replace- 
ment of gauge fields by matrices living on the links between lattice points. In this case, a 
U(l) covariant derivative will be substituted by 

Dj<f>{x) = j (e~ iLY ^ x U(x + xj) - </>(xf) . (A.2) 

In order to check whether our conclusions still hold using link variables, we first per- 
formed a series of simulations in the semilocal case to determine whether the evolution was 
compatible. Beginning from the same initial conditions, we evolved the system according 
to each discretization scheme. We observed that the simulations undergo extremely similar 
evolution on a pointwise basis, as seen in figure |4|, though there are modest differences 
with slightly longer strings on the lattice gauge calculations, occasionally leading to extra 
connections. 

To quantify this small difference and its impact on our statistical result, we performed a 
further 120 simulations using both discretizations. To save computing power we simulated 
our system in the semilocal case (sin 2 #w = 1)> discretized in 64 3 cubes for different values 
of (5. We used several different initial conditions, but for each one we let the system 
evolve using both the naive discretization of the (lagrangian) equations of motion and the 
link variable method for discretizing the (hamiltonian) equations of motion. The initial 
conditions were calculated using the method b) described in the text (section |3|), and as 
in the electroweak simulations performed in this work, we added an ad hoc damping term 
(77 = 0.5). As our interest lies mainly in the late-time behaviour of the dumbbell network, 
because the smaller cube size forces t < 64, we chose to calculate the total string length at 
t = 56 for both cases. 

Figure || shows the results. We computed the string lengths as explained in section 
footnote ||, and then added all string lengths for strings longer than five times their length. 
That final measure is represented in figure [B|, which shows clearly that the differences 
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Figure 5: Plot of string lengths in units of string widths (see text) for the naive discretization 
(squares) and the lattice link variable formalism (circles). The data are measured at time t = 56. 
The error bars are 1 a over 10 simulations. 

between the schemes are well within the uncertainties. Although performed for only 64 3 
cubes, this result should continue to hold for the larger cubes used to derive the main 
results in the body of this paper. The system simulated in this test is only the semilocal 
case, not the full electroweak case, but we are convinced that using the lattice link variable 
method will not alter significantly the results presented in this paper. 
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